UT-847 



End-point of the Electroweak Phase Transition using the 

auxiliary mass method 

Kenzo Ogure and Joe Sato* 
Institute for Cosmic Ray Research, University of Tokyo, Tanashi, Tokyo 188-8502, Japan 
* Department of Physics, School of Science, University of Tokyo, Tokyo 113, Japan 

Abstract 

We study the end-point of the Electroweak phase transition using the auxil- 
iary mass method. The end point is mu ~ 40 (GeV) in the case mj = (GeV) 
and strongly depends on the top quark mass. A first order phase transition 
disappears at ~ 160 (GeV). The renormalization effect of the top quark is 
significant. 

The Electroweak phase transition is one of the most important phase transitions at the 
early universe since it may account for the baryon number of the present universe [|I|. This 
phase transition was first investigated using the perturbation theory of the finite temperature 
field theory and predicts the first order phase transition from an effective potential [0|§. The 
perturbation theory, however, has difficulty due to an infrared divergence caused by light 
Bosons and cannot give reliable results in the case where the Higgs Boson mass m#, is 
comparable to or greater than the Weak Boson mass. Lattice Monte Carlo simulations, 
therefore, become the most powerful method and are still used to investigate details of the 
phase transition ||-§|. According to these results, the Electroweak phase transition is of the 
first order if m# is less than an end-point w,h,c ~ 70 (GeV). It turns to be of the second 
order just on the end-point. Beyond the end-point, we have no phase transition, which 
means any observable quantities do not have discontinuities. As far as we know, three other 



non-perturbative methods predict the existence of the end-point [pHTl|. The end-point is 
determined below 100 GeV by these three methods. 

The auxiliary-mass method is a new method to avoid the infrared divergence at a finite 
temperature T fT2}-|i"5fl. This method is based on a simple idea as follows. We first add 
a large auxiliary mass to light Bosons, which cause the infrared divergence, and calculate 
an effective potential at the finite temperature. Due to the auxiliary mass, the effective 
potential is reliable at any temperature. We next extrapolate this effective potential to 
the true mass by integrating an evolution equation, which we show later. We applied this 
method to the ^-invariant scalar model and the 0(iV)-invariant scalar model, and obtained 



satisfactory results ||T3HT5|]. 

We apply the method to the Standard Model and investigate the Electroweak phase 
transition in the present paper. We add an auxiliary mass M > T only to the Higgs Boson, 
which becomes very light owing to a cancellation between its negative tree mass and positive 
thermal mass for small field expectation values around the critical temperature. We notice 
that the infrared divergence from the Higgs Boson is always serious if the phase transition 
is of the second order or of the weakly first order [[J. In the standard model, transverse 
modes of the gauge fields also have small masses at small field expectation values since they 
do not have the thermal mass at one loop order. It is however, expected that they do have 
a thermal mass (~ g 2 T) at the two loop order [f7|P,|l6| . Here, g is a gauge coupling constant. 
If so, the loop expansion parameter || is < 1, even if the field expectation value is zero. 
Here, Mq is the mass of the gauge Boson, which is a sum of a zero-temperature mass and 
a thermal mass. We assume that this actually occurs and the infrared divergence from the 
gauge Bosons is not serious. Since this small thermal mass for the transverse modes will 
bring only a slight change to a one-loop effective potential, we use the one-loop effective 
potential without this small mass for the transverse modes. 

An effective potential is then calculated as follows in the Landau gauge @,[J, 
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In the above equations, A, #2, #1 and gy are coupling constants for the Higgs Boson, SU(2) 
gauge field, U(l) gauge field and top Yukawa respectively. The matrix T is orthogonal 
and diagonalizes the mass matrix for the Z Boson and photon at finite temperature. We 



renormalized the effective potential using the MS scheme with a renormalization scale fx. 
A zero-temperature contribution from the Higgs Boson is neglected since it is small in the 
mass region we consider. The ring diagrams are added only to the Weak Bosons and the 
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Z- Boson since the Higgs Bosons have auxiliary large mass and do not need the resummation. 

We then extrapolate this effective potential at the auxiliary mass squared M 2 to that of the 

true mass squared — v 2 using an evolution equation. Since we add the auxiliary mass only 

to the Higgs Boson, the evolution equation is same as that for 0(4)-invariant scalar model, 

which was constructed in[] 

dV 1 72 1 r°° „ k 2 1 
~(p H — n / dk- 



A non-perturbative effective potential free from the infrared divergence can be obtained by 
solving the evolution equation (0) with an initial condition Eq.(|I|) numerically. 

Before showing our numerical results, we relate the parameters z/ 2 , A, #2, Qi and gy to 
physical quantities at the zero-temperature |§, 
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Radiative corrections at the one-loop order are included in the equations for v 2 and A since 
they are large, especially in the case where the Higgs Boson mass is small. The effective 
potential Eq.([I]) does not depend on p, using A in Eq.(|3|) in this order. We fix the masses of 
the Weak Bosons and the Z-Boson as Mwo — 80 (GeV) and Mzo = 92 (GeV) below. 



1 We neglected the momentum dependence of a full self-energy in Ref. [14]. This corresponds to 
the local potential approximation of the systematic derivative expansion of the effective action. 
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We first investigate a SU(2) x U(l) gauge plus Higgs theory, corresponding to the case 
m t = 0. We show results obtained by setting M = T since similar results were obtained 
by setting M — -| and M = 2T as in the case of This is quite natural since the 

ristriction on M is M > T. The effective potentials at the critical temperature are shown 
in Fig.|l| for mn = 15, 30, 45 (GeV), respectively. The first order phase transition becomes 
weaker for smaller values of the Higgs mass and disappears finally They are compared to 
effective potentials obtained by the ring resumed perturbation theory at the one-loop order 
without the high temperature expansion in Fig.[| We find clearly that they are similar for 
smaller values of rnu and different for larger values of tuh- This is consistent with the fact 
that the ring resumed perturbation theory is reliable only for smaller values of the Higgs 
mass mn <C Myy ||. We plot a ratio of the critical field expectation values to the critical 
temperature, C /T C , as a function of m H in Fig|| This quantity indicates the strength of 
the first order phase transition and important in estimating the sphaleron rate, which plays 



a very important role in the Electroweak Baryogenesis JT7|,[T8| . The end-point is determined 
as rriH,c = 38 (GeV) from Fig.[5| This figure also shows that the results obtained by the 
auxiliary mass method and the perturbation theory is similar for smaller values of mn and 
different for larger values, m H > 30 (GeV). 
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FIG. 1. The effective potentials at the critical temperature for run =15, 30, 45 (GeV). The 
first order phase transition becomes weaker for smaller values of the Higgs mass and disappears 
finally. 
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FIG. 2. The effective potentials at the critical temperature obtained by the auxiliary mass 
method and the perturbation theory for ran =15, 30, 45 (GeV). They are similar for smaller 
values of run and different for larger values of run- 
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FIG. 3. The ratio of the critical field expectation values to the critical temperature, 4> C /T C . 
The results obtained by auxiliary mass method and the perturbation theory are similar for smaller 
values of mn and different for larger values, mn ^ 30 (GeV. 

We next investigate more realistic cases in which the top quark mass is finite. The same 
ratios are shown in Fig.|] for various values of m t . This figure shows that the strengths of 
the first order phase transition are almost same for m t 100 (GeV) and become weaker for 
m t ^ 100 (GeV) rapidly. The end-points are then shown in Fig.[5] as a function of m t . The 
graph labeled "1-loop" is obtained using Eq.(P and Eq.(^|), which take into account the zero- 
temperature radiative corrections from the top quark and gauge fields. The contribution 
from the top quark is much larger than that of the gauge fields. On the other hand, the graph 
labeled as "tree" is obtained without the zero-temperature radiative correction, omitting the 
contributions from foo and Jfo from Eq.(|l|) and leaving only the first terms of Eq.(|3]) for 
A and v 2 . They are not much different for smaller values of the top quark mass, m t < 100 
(GeV). Their behavior, however, differs drastically for larger values of the top quark mass, 
m t > 100 (GeV). Surprisingly, the end-point vanishes for m t > 160 (GeV) in the "1-loop" 
results though it increases in the "tree" results. These results tell us that fermionic degrees 
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of freedom play significant roles in the phase transition through the renormalization effects 
at the zero-temperature. We also conclude that there are no the first order phase transitions 
for m t = 175 (GeV), no matter how small the Higgs Boson mass. 
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FIG. 4. The ratio of the critical field expectation values to the critical temperature, cj) c /T c . The 
results obtained by the auxiliary mass method and the perturbation theory are similar for smaller 
values of mn and different for larger values, tuh <i 30 (GeV). 
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FIG. 5. The end-points as a function of mj. The graphs labeled as "1-loop" and "tree" are 
obtained with and without the zero-temperature radiative corrections from the top quark and 
gauge fields, respectively. The end-point vanishes for m t ^ 160 (GeV) in the "1-loop" result. 

In the present paper, we have calculated the effective potentials of the standard model 
using the auxiliary mass method at a finite temperature. We first investigated a SU(2) x U(l) 
gauge plus Higgs theory, corresponding to the case m t = 0. The phase transition was of 
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the first order and similar to the results obtained by the perturbation theory for smaller 
ran ~ 15 (GeV). The phase transition became weaker for larger Mn ~ 30 (GeV) and finally 
disappeared in contrast to the results from perturbation theory. We found that the end- 
point is at rriH,c — 38 (GeV) in this case. This is consistent with the results of the Lattice 
Monte Carlo simulation |5|-|| and the other non-perturbative methods PHTT| qualitatively. 
The value of the end-point, however, was smaller than those by these methods. This may 
be caused by the approximations, used to construct the evolution equation @, or used in 
the other papers. The two loop effect from the gauge fields may shift our results due to 
slow convergence of the perturbation theory. We next investigated the more realistic case 
in which the top quark mass is finite. We found that the end-point was strongly dependent 
on m t and disappeared for m t > 160 (GeV). The renormalization effects from the top quark 
were significant. Lattice Monte Carlo simulations, however, do not follow this behavior [|J. 
We think of two possible reasons: (1) Since our results differ from that of the Lattice Monte 
Carlo simulation by factors of 2 in a SU{2) x U(l) gauge plus Higgs theory quantitatively, 
the similar behavior may be found at a larger top quark mass in the Lattice Monte Carlo 
simulation. (2) Since the one-loop correction to the effective potential at the zero-temperature 
is significant, the 3D effective theory, which has no Fermionic degrees of freedom, may not 
reflect the effect appropriately. 

Finally, the strongly first order phase transition necessary for the Electroweak Baryoge- 
nesis was not found in the Standard Model. We will apply this method to extensions of the 
Standard Model. 
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